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ABSTRACT 

We formulate a general, steady-state model for the torque on a magnetized star 
from a surrounding accretion disc. For the first time, we include the opening of dipolar 
magnetic field lines due to the differential rotation between the star and disc, so the 
magnetic topology then depends on the strength of the magnetic coupling to the 
disc. This coupling is determined by the effective slip rate of magnetic field lines that 
penetrate the diffusive disc. Stronger coupling (i.e., lower slip rate) leads to a more 
open topology and thus to a weaker magnetic torque on the star from the disc. In the 
expected strong coupling regime, we find that the spin-down torque on the star is more 
than an order of magnitude smaller than calculated by previous models. We also use 
our general approach to examine the equilibrium ('disc-locked') state, in which the net 
torque on the star is zero. In this state, we show that the stellar spin rate is roughly 
an order of magnitude faster than predicted by previous models. This challenges the 
idea that slowly-rotating, accreting protostars are disc locked. Furthermore, when the 
field is sufficiently open (e.g., for mass accretion rates ^ 5 x 10~^Mq yr~^, for typical 
accreting protostars), the star will receive no magnetic spin-down torque from the disc 
at all. We therefore conclude that protostars must experience a spin-down torque from 
a source that has not yet been considered in the star-disc torque models — possibly from 
a stellar wind along the open field lines. 
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1 INTRODUCTION 

Accretion discs are responsible for some of the most en- 
ergetic and spectacular phenomena in many classes of as- 
trophysical objects, including protostars, white dwarfs (cat- 
aclysmic variables and intermediate polars), neutron stars 
(binary X-ray pulsars), and black holes (both stellar mass 
X-ray transients and supermassive active galactic nuclei). 
Gravitational potential energy liberated by the accretion 
process gives rise to exceptional luminosity excesses and can 
drive powerful jets and outflows. Accretion onto the central 
object can occur only as quickly as angular momentum can 
be transported away from the system. Furthermore, the ac- 
cretion of disc material, which has high specific angular mo- 
mentum, spins up the central object, if the object rotates at 
less than the break-up rate. It is therefore surprising that 
the central objects (hereafter 'stars') are often observed to 
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spin far below their break-up rates, in spite of long-lived 
accretion. Why does this happen? 

There is good evidence that accretion onto magnetized 
stars occurs along closed magnetospheric field lines that con- 
nect the star to the inner edge of the disc. For instance, the- 
oretical models of this sort have been successful in explain- 



'Havashi et al.l Il996l: iGoodson et al,. 
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200l|), intermediate polars (e.g., 
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pulsars (e.e..|Joss & Rappapord 
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evidence that the stars are m agnetized, namely for accret- 
ing protostars JJohns-Krull et aL..1999.'l , in termediate polars 
(Pii rolaet alj[^93h . and X-rav pulsars (IMakishima et al] 

Magnetic fields can also be effective at transferring 
angular momentum away from the star, possibly explain- 
ing the observed rotation rates. Torques on the star that 
are exerted by magnetic field lines anchored to the star 
and that are also connected to the disc have been cal- 
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culated by several authors (e.^.. [c hosh fc LambI Il979l: 
Cameron&^an^belj[ 119931: iLovelace et all Il99a IWand 
95': 'Yi"l995': 'Armit aia-e fc Clarkel Il99fil . hereafter AC96; 
Rappaport et a l. 2004). Under certain circumstances, this 



i iply explamm K tne slow spin ol some iU^^^Uiyger^_^y^J]Newm^iet_a!J litJij^ : IL/Ovelace et 
lKonig]|[l99ll . hereafter K91). In this ll995l:lBardou fc Hevvaert Jl996ll . 'IVpicallv. this inflation 



torque can counteract the angular momentum deposited by 
accretion, leading to a net spin-down of the star (possi- 
bly e xplaining spin-down e pisodes observed in X-ray pul- 
sars; iGhosh fc Laiiibl 119781: ILovelace et aLlll995D or giving 
rise to an equilibrium state, in which the net torque on 
the star is zero (poss ibly explainin g the slow spin of some 
accreting protostars; 
equilibrium state, the spin rate of the central object de- 
pends on the accretion rate in the disc, and so a system 
is then considered to be 'disc locked.' Since these mod- 
els for the magnetic star-disc interaction show that accret- 
ing stars can spin more slowly than the break-up rate, 
there is a general perception that the presence of an ac- 
cretion disc in any system leads to slow rotation rates. 
This idea of disc locking has been applied to a variety 
of problems. As an example, in systems where the mo- 
ment of inertia of the star is changing (e.g., during contrac- 
tion) , some authors have assumed that disc- locking keeps the 
st ar at a constant spin rate (e.g., as applied to protostars 



after UKL). They have shown that, as the differential twist 
angle between the star and disc A$ monotonically increases, 
the torque exerted by field lines first reaches a maximum 
value, then decreases. This occurs because azimuthal twist- 
ing of the dipole field lines generates an azimuthal com- 
ponent to the field, and the magnetic pressure associated 
with this component acts to inflate the held, which then bal- 
loons outward at an angle of ~ 60° from the rotation a xis, 
causally disconnecting the star and disc (se e also lAlvlil98a: 
Alv fc Kuiiperslll990t iNewman etlo] Il992l: iLovelace etlo] 



St ar at a constant spin rate (e.g.. as applied to protostars 
bvlBouvier et aLlliggTl: ISills et alil2000HBarnes et al.ll200ll: 
iTinker et alJ2002l : lRebull^et alJ2004l: and sugges t ed for stel- 
lar collision products by 'L eonard fc Livid [l995l : ISills et all 
I2OOI: De Marco et al. 2004) . 



There is a nagging problem with this physical picture, 
however, because the magneti c torque calculations d iscussed 
above (with the exception of ILovelace et al']|l99 5') assume 
that the stellar magnetic field remains largely closed and 
that field lines connect to a large portion of the disc^. This 
assumption is questionable because closed magnetic struc- 
tures tend to open when enough energy is added to them, 
and these systems possess a natural source of energy in the 
form of gravitational potential energy that is released dur- 
ing disc accretion. This energy release can drive outflows and 
twist fleld lines, thereby adding energy to the magnetic held. 
Thus, the general surplus of energy in accreting systems sug- 
gests that associated magnetic flelds should be dominated 
by open, rather than closed, topologies. How are low spin 
rates achieved in this case? 

In this paper, we generalize the star-disc interac- 
tion model to include the effect of varying field topol- 
ogy (i.e., connectedness). We consider the mechanical en- 
ergy that is added to the field via differential rotation be- 
tween the star and disc as the only mechanism responsi- 
ble for opening the field (though our formulation is eas- 
ily adaptable for other mechanisms). The time-dependent 
behavior of a dipole stellar field attached to a rotating, 
conducting disc has been stud ied, using an analytic ap- 
proach, by several a uthors Ce. g.. iLvn den-Bell fc Boil ylll994l : 
lAgapitou fc Papaioizou ,200ol : juzdeng^^F£^n2002ar here- 



Note that the X-wind model o f lShu et al.l il994l. and subsequent 
works) is unique among the star-disc interaction theory in the 
literature, since it assumes that a system will always accrete very 
near its disc-locked state. The magnetic field geometry employed 
by the X-wind model is also unique and was designed, in part, to 
avoid the problem of field opening due to differential twisting, as 
considered in this paper. Therefore, much of our discussion does 
not apply to the X-wind. 



opening of the field occurs when a critical differential rota- 
tion angle of A"I>c ~ ti" has been reached, and the amount of 
flux that opens depends on the strength of the magnetic cou- 
pling of the field to the disc (UKL). This analytic work on 
the field opening has been corroborated by time-dependent, 
numerical magnetohydrodynamic simulations of the stellar 
dipol e -disc i nteraction fa ayashi et al. 1996; Goodsonet^lJ 
| l997l I199A iMiller fc Stona il997t iKato et ab i2001l l200i 
iMatt et alJl2002l : lRomanova et allbOOal : iKiiker et alJl20m 

Our primary goal in this paper is to determine the ef- 
fect of the topology of the magnetic field on the torques in 
the st eady-state, star-disc in teraction model. In a previous 
paper jMatt fc Pudrit J2004^ ■ we gave a brief outline of this 
theory and showed that a more open (i.e., less connected) 
field topology results in a spin-down torque on the star that 
is less than for the closed field assumption. Consequently, 
the equilibrium state (with a net zero torque) features a 
faster spin than predicted by previous models, which calls 
into question the general belief that accretion discs neces- 
sarily lead to slow rotation. The present paper contains a 
more detailed presentation of the theory and our assump- 
tions, and we consider all possible spin states of the system 
(not just the equilibrium state). We also extend our analy- 
sis to show that there are at least three different modes in 
which a magnetic star-disc system can operate. Our anal- 
ysis is applicable to all classes of magnetized objects that 
accrete from Keplerian discs. However, since an abundance 
of observational data exists for accreting protostars, in par- 
ticular for classical T Tauri stars (CTTS's), we adopt a set 
of fiducial parameters that are appropriate for these systems 
and discuss various aspects of the model in this context. 

Section 121 contains a formulation of the general model. 
The special case of a disc locked system is the topic of section 
|H] The final section (S^J contains a summary of our results 
and includes a list of problems with using the disc-locking 
scenario to explain CTTS spins, plus a discussion of three 
possible configurations of the general system. 



2 STAR-DISC INTERACTION MODEL 

Magnetic, star-disc interaction models in the literature differ 
in their various assumptions, adopted parameter values, and 
in the introduction of 'fudge factors, ' but they are quite 
similar on the whole (for a review, see lUzdenskvlliooil . We 
formulate a general model that builds upon this previous 
work (mostly following AC96), by including the effect of 
varying magnetic field topology, via the introduction of the 
physical parameters /? and 7c (defined below). 

According to the usual model assumptions, a rotating 
star is surrounded by a thin, Keplerian accretion disc. The 
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angular momentum vector of the disc is aligned with that of 
the star, which rotates as a solid body and at a rate that is 
some fraction of break-up speed, defined by 

where 57*, i?*, and A/* are the angular rotation rate, radius, 
and mass of the star, respectively. Note that / is always 
within the range from zero to one^. The disc rotates at a 
different angular rate than the star at all radii, except at 
the singular corotation radius given by 

Rco = r^^^'R*. (2) 

For r < _Rco, the disc rotates faster than the star, while for 
r > iico, the angular rotation rate of the star is greater than 
that of the disc. 

The disc is assumed to be in a steady-state wherein the 
mass accretion rate Ma is constant in time and at all radii. 
In a real disc from which winds are launched. Ma may have a 
weak radial dependence, but we assume this has a negligible 
effect on the model. A rot at ion- axis- aligned dipole magnetic 
field, anchored into the stellar surface, also connects to the 
disc. The field is strong enough to truncate the disc at some 
inner location Rt from where disc material is subsequently 
channeled along magnetic field lines as it accretes onto the 
star. In general, the disc may have its own magnetic field 
(either generated in a disc dynamo or carried in by the disc 
from larger scales). We do not specifically include this field 
in the model, though it may be responsible for angular mo- 
mentum transport in the disc (providing Ma) and may also 
aid in the connection of the stellar field to the disc. Within 
the disc, the kinetic energy of the gas is much greater than 
the magnetic energy of the stellar field, but the region above 
the star and disc (the corona) is filled with low density ma- 
terial, and so the corona is magnetically dominated. 

In this configuration, the magnetic field connects the 
star and disc by conveying torques between the two. Torques 
are conveyed on an Alfven wave crossing time, which is much 
shorter than the Keplerian orbital time. Everywhere that the 
magnetic field connects the star to the disc, exept at Rco, the 
magnetic field is twisted azimuthally by differential rotation 
between the two. Inside Rco the field is twisted such that 
field lines 'lead' the stellar rotation, so torques from field 
lines threading the region r < Rco act to spin up the star 
(and spin down the disc). Conversely, torques from field lines 
threading r > _Rco act to spin the star down (and spin the 
disc up). The accretion of disc matter onto the star also 
deposits angular momentum onto the star. 

In order for disc material to accrete, Rt must be less 
than Rco so that accreting material loses angular momen- 
tum to the star as it falls inward. In order for the star with 
/ < 1 to feel any spin down torques from the disc, the stel- 
lar field must connect to the disc beyond Rco- Under this 
condition, the field that connects outside Rco transfers an- 
gular momentum from the star to the disc. To maintain 
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Figure 1. Magnetic star-disc interaction. The stellar field con- 
nects to a region of the disc, from _Rt to -Rout, reaching beyond 
-Rco in this case. The stellar field dominates the accretion flow 
onto the star (arrow). 



a steady accretion rate, the disc must then transport this 
excess angular momentum outward, resul ti ng in an altered 
disc structure jSu nvaev fc Shakurall977al lbt lSpruit fc TaamI 
ll993l : lRappaport et al We will define -Rout as the out- 

ermost radial extent of the magnetic connection, and the 
usual assumption is that -Rout ^ ^^co . Figure Q illustrates 
the basic picture, and shows the locations of -Rt, Rco, and 
7?out for a possible configuration of the star-disc interaction 
model. 

The assumption of a dipole field refers to the poloidal 
component of the field, -Bp = {Bf. + BIY^'^ , where Br and 
Bz are the cylindrical r and z components of the magnetic 
field, and the closed field only exists in the region interior 
to 7?out. The dipole field is used for simplicity and because 
it has the weakest radial dependence [Bz oc along the 
equator) of any natural magnetic multipole. In reality, the 
twisting of dipole field lines alters the poloidal field, but 
this perturbation should be slight in the region where the 
magnetic field remains closed (as justified by the wo rk cited 
in 31] e.g., UKL, and see also lLivio fc PringldlTogl) . 

For discussion throughout this paper, it is often instruc- 
tive to use physical units, especially for comparison with ob- 
servations. For this purpose, we adopt a set of observation- 
ally determined 'fid ucial parameters' that are ap propriate 
for CTTS's (e.g., see I Johns-KruU fc Gaffor'dll2002l) : 

Ma = 5 X 10"** Mq yr"\ 
M, = 0.5 Mq, 

i?* 2 Rq , 

B, = 2 X 10^ G, 

where -B* is the stellar magnetic field strength at the equa- 
tor"^. However, our formulation of the problem is applicable 
to any magnetic star-disc system. 



Disc accretion soluti ons do exist for / > 1 am fc NaravanI 

Il99lt IPaczvnskilll99ll) . in which the star is actually spun down 
by accretion toward / = 1, even without any magnetic torques. 
However, we only consider cases with / < 1 in this paper, since 
this characterizes the spin of observed protostellar systems. 



^ These fiducial parameters arc slightly different than used for 
figures 2 and 3 of .Matt & Pudritz 1 2004) . who considered the 
specific case of the CTTS BP Tau. 
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2.1 Twisting and slipping of magnetic field lines 



The torque exerted by magnetic field lines threading an an- 
nulus of the disc of radial width dr is given by (e.g., AC96) 

,2 



where 



7 = B^/S, 



(3) 



(4) 



Here, ^ is the dipole moment and is the azimuthal com- 
ponent of the magnetic field. The radial component, Br, is 
assumed to be negligible within the disc, the torque has been 
vertically integrated through the disc, and B^ refers to the 
value at the disc surface. Here, and throughout this paper, 
we choose the sign of the torque to be relative to the star 
such that a positive torque spins the star up, and conse- 
quently spins the disc down (and vice versa for a negative 
torque) . 

The differential magnetic torque of equation |3| depends 
not only on ^, but also 7, which is the 'twist*,' or pitch angle, 
of the field. The total (integrated) magnetic torque will also 
depend on the size and radial location of the magnetically 
connected region in the disc. While /i is a constant parameter 
of the system, the radial dependence of 7 depends on the 
physical coupling of the magnetic field to the disc. 

In general, the coupling is not perfect. Magnetic forces 
act to resist the twisting of the field, and so the field will 'slip' 
backward through the disc at some rate Ud proportional to 7. 
The exact slipping rate depends upon which physical mech- 
anism is at work. In the literat ure, there ar e generally three 
mechanisms discussed (e.g., see lWanJl995l) : 1) magnetic re- 
connection in the disc, 2) reconnection outside the disc, and 
3) turbulent diffusion of the magnetic field through the disc. 
We adopt the latter mechanism, but we further discuss the 
other two, below. 

The magnetic fi eld slips azimuthally at a speed (e.g., 
iLovelace et al]ll995l: UKL) 



* This should not be confused with the 'twist angle' of the foot- 
points of the field, A^, as discussed by UKL, though A<I> and 7 
are intimately related. 




Vd = ^7 = /3i;k7, (5) 

where ryt is the turbulent magnetic diffusivity, h is the local 
disc scale height, Vk is the Keplerian orbital speed, and we /3 ~ 10 
have introduced the dimensionless 'diffusion parameter' /3 = 
?7t(/ii'k)~^- The variable /3 simply parametrizes the coupling 
of the stellar magnetic field to the disc such that /3 ^ 1 
corresponds to weak coupling, and /3 ^ 1 to strong coupling. 

Generally, /3 is a scale factor that compares to Wk, 
and we have chosen this generic formulation so that the sys- 
tem behavior is largely independent of any particular disc 
model — as long as the disc obeys Keplerian rotation and 
provides a steady accretion r ate. However, if we tempor arily 
consider a standard a-disc dShakura fc Sunvaevlll9'73l . we 
may rewrite our diffusivity parameter /3 in a more physi- 
cally revealing way: 

where a has its usual meaning and Pt is the turbulent 



Figure 2. The magnetic twist (7 = B^/B^) times the diffusion 
parameter /3, as a function of radius along the surface of the disc 
(thick, solid line). The dotted lines indicate critical radii for a 
field-opening twist of 7c = ±0.5//3 (chosen for illustrative pur- 
poses). The locations of Ri^, Rco, and Rout are where 7 equals 
—7c, 0, and 7c, respectively. 



Prandtl number, equal to the turbulent viscosity divided by 
rjt . The disc turbulence is lik ely to be driven by the m agneto- 
rotational instability (MRI: lBalbus fc Hawlev|ll99 j) . which 
follows the general behavior of an a-disc. Since both a and 
h/r typically have weak radial dependences, and the value 
of Pt is unknown, we assume that (3 is constant in the region 
of the disc connected to the stellar field. 

The value of (3 is not well constrained (AC96 used /3 — 
1), but extreme a-disc parameters give an upper limit of 
/3 < 1. For a more reasonable estimate, note that a thin disc 
usually mea ns /i/r<0.1j and a is typically in the range of 
0.001 to 0.1 JSano et alil2004^ . So, assuming Pt is of order 
unity, f3 < 10"'^. We get a similar estimate using equation |K| 
and reasonable guesses for CTTS disc parameters. 



. 10" cm2 



h 

Rq 



«k 



100 km s- 



(7) 



However, given the uncertainties and possible variation 
among different systems, and to assess the effect of the cou- 
pling of the field to the disc, we retain /3 as a free parameter. 

In the disc connected region, if Vd is anywhere less 
(greater) than the local differential rotation speed between 
the star and disc, 7 will increase (decrease) on an orbital 
time-scale. Thus the magnetic field will quickly achieve a 
steady-state configuration in which Vd equals the local dif- 
ferential rotation rate (e.g., UKL), which gives 

7 = /3-'[(r/Pco)'/'-l]. (8) 

The solid line in Figure|2|shows the quantity [3"f as a function 
of radius (normalized to Pco) along the surface of the disc. 
The magnetic twist is zero at Pco, and is oppositely directed 
on either side of Pco. Also, at a given radius, the twist will 
be larger for smaller values of /3 (and vice versa). 

We have assumed that the field coupling is determined 
by turbulent diffusion, and when f3 is constant, we find that 
7 oc r^'^ (eq.|H| for r 2> Pco). Other coupling mechanisms (as 
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discussed ab ove) result in a different radial dependence of 7. 
For example. lWand (Il995l) showed that if the twist is limited 
by reconnection in the disc, 7 oc r^'^ (for r ^ Rco), while 
for reconnection outside the disc, 7 appro aches a constant 
value (for r ^ Rco)- On the other hand. IlIvIo fc Pringld 
il992h and AC96, assumed 7 was limited by reconnection in 
the stellar corona, and they used the same formulation as 
equation |H| (with P — 1). In any case, note that the radial 
dependence of the differential magnetic torque is dominated 
by the falloff of the dipole magnetic field (which results in 
the r~'^ dependence of eq.|^. Therefore, the choice of mag- 
netic coupling mechansim will not much inffuence our results 
(AC96). Similarly, a small radial dependence of P (which we 
take as constant) will not introduce a large error. 



2.1.1 Maximum twist for dipole field 

As discussed in section several authors have shown that 
dipole field lines will transition from a closed to an open 
topology when a critical differential rotation angle of A$c ~ 
TT has been reached. This corresponds to a critical field twist 
of 7c ~ 1. Since the twisting of field lines does not signifi- 
cantly alter the poloidal configuration of the field lines that 
remain closed, and as an approximation, we will assume that 
the opening of field lines is only important for the determi- 
nation of the size of the connected region (i.e., to determine 
Rout). Thus, we include the effect of field line opening in 
the steady-state torque theory in the following manner: we 
will use equation |H| only where 7 < 7c and assume that the 
field will be open everywhere else (a similar approach was 
used by iLovelace et alJ Il995l and justified by the work of 
UKL). In other words, wherever equation |H| predicts 7 > 7c, 
the magnetic connection is assumed to be severed, so the 
star and disc are causally diconnected, such that no torques 
can be conveyed between the two. The size of the connected 
region in a Keplerian disc is then limited to a finite radial 
extent near Rco, where the differential rotation between the 
star and disc is the smallest. 

Will field lines, once opened by differential rotation, re- 
main open? It has been suggested that such field lines could 
reconnect in the current s heet formed during the opening 
jAlv fc Kuiipersl lT990: Uzd enskv erai]l2002lJl . In order for 
this to be important, the time-scale for reconnection should 
be comparable to or shorter than that for field line open- 
ing. I t is not clear whether this is the case in these sys- 
tems dUzdenskv et alj|2002bt iMatt et alJl2002l) . but even ff 
it is, there are other considerations. First, due to the topol- 
ogy of the field, reconnection must initially occur between 
open field lines at the smallest radii (connecting to the low- 
est latitude on the star). It is possible that, if reconnec- 
tion does occur, only a small amount of fiux will be able 
to reconnect before this newly co nnected field again b egins 
to open (as in the simulat ions of iGoodson et aOll999l and 
see lUzdenskv et~ai]l2002b^ . In this case, the size of the con- 
nected region (in a time-averaged sense) will be only slightly 
larger than if the reconnection were never to occur. Second, 
the configuration of the opened field is favorable to launch 
magnetocentrifugal flows from the disc iBlandford fc Pavnd 
Il982r . We ignore such outflows in our model, but in a real 
system, they could help to maintain an open magnetic field 
configuration. Thus we conclude that, once the field has 



opened, reconnection along the current sheet is unlikely to 
significantly affect the size of the connected region. 

2.1.2 Maximal spin-down torque for maximal twist 

It is instructive to look at the maximum possible spin-down 
torque in this system. Regardless of any disc model or any 
magnetic coupling physics, the largest possible magnetic 
torque on a star that connects to a disc via a dipole magnetic 
field occurs when the field is maximally twisted (7 — 7c) at 
all radii along the surface of the disc. Spin-down torques 
on the star only occur along field lines threading the disc 
outside iico. Also, the accretion of mass from a Keplerian 
disc always adds angular momentum to the star. Therefore, 
the largest net spin-down torque on the star occurs when the 
disc is truncated exactly at the corotation radius {Rt = -Rco), 
the field threads the disc to -Rout —> 00, and the disc does 
not accrete (-Ma = 0). One then integrates equation |3] from 
7?co to CO to get 

max To t^ /q\ 

'twist 2 ^3 ■ V / 

This is the absolute maximum spin-down torque that a star 
can undergo from a disc that exists in the equatorial plane 
and to which the star is connected via a dipole magnetic 
field. It is even independent of the rotation profile of the 
disc, except that angular rotation rate of the disc is slower 
than the star outside some radius Rco- It is also independent 
of the angular momentum transport mechanism within the 
disc. 

To achieve this maximal torque requires that a) the 
field twist has no radial dependence and b) the twist is very 
nearly equal to the maximum allowed value of 7c. If the 
coupling of the field to a Keplerian disc is determined by 
turbulent diffusion, a constant 7 can only be achieved in 
the unlikely event that rjt decreases with radius to exactly 
counteract the increase in differential rotation rate. Alter- 
natively, reconnectio n in the stellar or disc corona may also 
lead to a constant 7 fAlv fc Kuiiperslll990l:IWandll995l but 
see discussion in i|2.1.H . However, in either case, it is not 
clear why the value of the constant twist would necessar- 
ily be near the maximal value 7c (instead of, e.g., 0.l7c). 
Though this torque may not be very realistic, it is similar 
in strength to the spin-down torque used in previous models 
(e.g., it is the equivalent to the solution of AC96, for 7c = 1). 

2.1.3 Determination of Rm and Rout 

To derive a more realistic magnetic torque, we adopt 
equation El for the radial dependence of 7. Following 
iLovelace et al.l il995l) . we assume that this equation is only 
valid where I7I < 7c, and that the field is open everywhere 
else. Thus, equation |H| predicts that the outer radius of the 
magnetically connected region in the disc is 

-Rout = (1 + 13-fc)^^'' Rco. (10) 

There is a corresponding location insided -Rco at which the 
twist formally exceeds the critical value, given by 

Riu = {l- p-ycf'Rco. (11) 

The dotted lines in Figure |2| indicate these radii for /37c — 
0.5, in which case -Rin ~ 0.63-Rco and -Rout ~ 1.31 -Rco. It 
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is evident that the field topology is a function of /37c such 
that more diffusion in the disc allows for a larger connected 
region. Note that, if /37c > 1, the field can remain connected 
to the disc at any r < 7?co (since i?in is then not defined). 

The typical assumption of a closed magnetic topology, 
corresponding to i?out oo (e.g., AC96; |Yj[1993), is equiv- 
alent to 7c — > oo — the field is allowed to twist to arbitrarily 
large values without opening. In order to consider the ef- 
fect of varying topology (i.e., where the field is open beyond 
some finite -Rout), we adopt a value of 7c = 1 (as justified 
by, e.g., UKL). However, we will retain 7c as a parameter in 
all of our formulae for a comparison between the two cases 
(7c — > oo and 7c — 1) and so that different values of 7c may 
be considered by the reader. The combined parameter Pjc 
appears throughout our formulation. We generally think of 
this parameter in two ways. First, when 7c = oo, the stellar 
field is closed and connects to the entire disc, and this rep- 
resents the 'standard' star-disc interaction model. Second, 
for the more realistic case of 7c = 1, the field topology is 
partially open, and -Rout then depends on (3. 

2.2 Three possible states of the system 

The inner edge of the disc is delimited by -Rt (discussed 
in tl2.3.1^ . and so there are two possible magnetic configu- 
rations of an accreting system, depending on the location 
of Rt relative to -Rin. First, if -Rt < -Rin, the stellar field 
will be largely open, and equation |H1 is not valid anywhere. 
The outer radius of the magnetically connected region, -Rout 
(which is not then determined by eq. 110^ will be near the 
inner edge of the disc (-Rout ~ -Rt). We will refer to this sit- 
uation as 'state 1' of the system. In state 1, the star receives 
no spin-down torques from the disc. 

In 'state 2,' -Rin < Rt < Rco and the star is magnetically 
connected to the disc from -Rt to -Rout. State 2 represents 
the typical configuration considered in many models and 
was discussed at the beginning of this section. Also, systems 
near their disc-locked state (f|2Jl are always in state 2. 

Finally, there exists a third possible, non-accreting 
state, 'state 3,' that occurs when the disc is overpowered 
by the magnetic field (e.g., for fast rotation, large /i, or 
small - Ma) and the disc bec omes truncated outside -Rco (e.g., 
[lUarion ov fc SunvaevllQTsh . In state 3, there are no positive 
(spin-up) to rques on the star, so it ca n never be in spin 
equilibrium JSunvaev fc ShakurallToTTbl) . 

Figure |3] illustrates the basic magnetic configuration of 
each state. One can think of this Figure, for example, as a se- 
quence (from top to bottom) of decreasing A/a (or increasing 
n) . A steadily decreasing -Ma may represent an evolutionary 
sequence (e.g.) for protostars as one goes from class sources 
to weak- lined T Tauri stars (class 3) . As -Ma decreases from 
a system in state 1, the disc truncation radius (which is at 
the inner edge of the disc in the Fig.) moves outward and 
eventually crosses the location of -Rin (entering state 2) and 
then -Rco (state 3). 

The conditions for which a system tr ansitions from an 
accreting state to state 3 is unknown (see lRappaport et alJ 
I2OO4I and t|4.2.3ll . In the current work, we do not consider 
state 3, other than to note that it occurs somewhere below 
the lower -Ma limit of accreting systems. Instead, we focus 
most of our attention on state 2 and discuss state 1, where 
appropriate. 




Figure 3. Three possible configurations of the magnetic star- 
disc interaction. There are two possible accreting states: either the 
stellar field connects only to the inner edge of the disc (top panel) 
or it connects to an extend region, reaching beyond -Rco (middle 
panel). In either case, the stellar field dominates the accretion 
flow (arrows) onto the star. Under certain conditions (e.g., low 
-Ma), accretion onto the central star will cease, defining the third, 
non-accreting state (bottom panel). 



2.3 Torques between the star and disc 

A combination of equations |3 and |H| gives the full radial 
dependence of the differential magnetic torque in the system. 



dr (5 Rt \Rco) \Rco) 



(12) 



where, for convenience, we have used equation |21 to express 
the radius in units of -Rco- Furthermore, the angular momen- 
tum carried by accreting material through each annulus of a 
Keplerian disc (of width dr and ve rtically integrated) equals 
dra = Q.^M^{GK'h/r)-^/^dr (e.g., |ciar£lertl][l99|). This 
can be combined with equation|5|to give the differential ac- 
cretion torque, as a function of r/Rco, 



dr, 



GM, 



1/2 



(rco) 



-1/2 



(13) 



The spin of accreting stars 7 



in 

c 



-a 





^ ^ : 




/ ^ ^internal : 




/ ~~ ^ ^ : 

-/^^. accretion ~~ - ^ _ ; 

'\ 


: Rt, 


\ magnetic^ ! 


: 1 
: 1 
:...(. 





0.5 1.0 1.5 2.0 2.5 3.0 



Figure 4. Differential torques in the fiducial CTTS system (see 
discussion above i|2.1l for /3 = 1 and 7c — > oo, where the field 
is assumed to remain connected over the entire disc. The trun- 
cation radius, i?t, is where dr^a = (ifa and drj = 0, indicated by 
the vertical dotted line (i?t ~ 0.91-Rco)- The system is shown in 
its equilibrium state, where the net torque on the star is zero, 
requiring a stellar spin period of 6.0 days (iJco ~ 5.5-R, ) 



The assumption that is constant at all radii in the 
disc, requires that the net angular momentum transported 
away from each annulus in the disc equals dra. The disc must 
therefore be structured in such a way that the differential 
torques internal to the disc, dri, satisfy 



dri = dra — drn 



(14) 



These internal torques could result fro m angular momentum 
trans port vi a (e.g.) turbulent viscosit y dShakura fc SunvaevI 
11973^ . MRI (iBalbus fc Hawlev^ll991^ . or disc winds (see re- 
view bv lKonid fc Pudritzl2000ir If one assumes a particular 
angular momentum transport mechanism in the disc, the so- 
lution to equation 1141 determines the stru cture of the disc. 
As an example, for the case of a viscosity, iRappaport et all 
l|2004fl showed that the disc can respond to external mag- 
netic torques by increasing its surface density in order to 
transport the additional angular momentum outward. A de- 
tailed treatment of the disc adjustment is not necessary here, 
since we are presently concerned with torques on the star, 
and we simply assume that the disc structures itself such 
that eg nation 1 141 is valid. 

Figure 2] shows the differential torques (dr/dr) as a 
function of r/Rco for the adopted fiducial parameters. The 
solid, dcish-dotted, and dashed lines in Figure 01 represent 
the differential torques from equations 1121 1131 and 1141 re- 
spectively. The system shown has (3 — 1 and 7c — > oo, so 
that the field is connected to the entire surface of the disc, 
and so that the Figure represents the closed topology of sev- 
eral models in the literature (e.g., AC96). 

The differential magnetic torque (drm, solid line in 
Fig. 2J is strongest near the star, where the dipole field is 
strongest, and it acts to spin up the star (and thus spins 
down the disc) for r < _Rco • At Rco, drm goes to zero, since 
the field is not twisted (7 = 0) there. Outside Rco, the mag- 
netic torque becomes stronger again, as the twist increases. 



though now acting to spin the star down (and the disc 
up). The dipole field strength falls off faster with distance 
{Bz oc r~^) than the magnetic twist increases (7 oc r^^^), so 
drm has a minimum value at r/Rco ~ 1.37 and then goes to 
zero as r ^ 00. 

The dashed line in the Figure (dri) gives us some infor- 
mation about the disc structure. In this case, the structure 
will be significantly different than for a case with drm = 0. 
For a very large drm outside iico, the assumption that the 
disc can counter act it (via an increase in dri) must eventu- 
ally break down, and the system would then be in state 3 

(ail- 



2.3.1 Truncation radius 

Inside the corotation radius, the stellar magnetic torque acts 
to extract angular momentum from the disc, further en- 
abling accretion (not hindering accretion, as for r > i?co). 
The differential magnetic torque increases rapidly as one 
moves toward the star, and at some point, drm = dra. There, 
the external magnetic torque alone is capable of maintain- 
ing Ma- Consequently, dri goes to zero at the same radius 
and formally becomes negative for smaller r (see Fig. UJ. 
Negative dri is unrealistic, however, since it would require 
angular momentum transfer inward through the disc, from 
slower spinning material to faster spinning material, so the 
Keplerian disc does not exist where where dri < 0. Thus, 
the disc truncation radius^, Rt, is where dn = 0. 

At _Rt, the stellar magnetic field will quickly spin down 
the disc material, forcing it into corotation with the star. 
Sub-Keplerian rotation leads to a free-fall of disc material 
onto the surface of t he star in a 'funnel fi ow' along magnetic 
field lines (e.g., K91: lRomanova et al.l2'o02) . Whether or not 
the funnel flow originates exactly f rom Rt or from insid e 
that radius is subject to debate (e.g.. lAlv fc Kuiiperj|l990ll . 
However, for the present discussion of angular momentum 
transport, the most important thing is that Rt defines the 
location where the stellar magnetic field completely domi- 
nates over the disc internal stresses, and so all of the angular 
momentum of disc material at Rt will end up on the star. 

By setting eauation ll2l eaual to ll3l we derive a relation- 
ship defining Rt , 



V Rco J V Rco J 



P_ r-7/3 

I J 5 



where 

= 2/i^Ma"'(GAf. 



1/2^-7/2 



(15) 



(16) 



is a dimensionless parameter relating the strength of the 
magnetic field to the strength of accretion. This formula was 
also derived by Yi ( 1995), but with different disc parameters 
in place of our (3. For any given (3, f, and ip, there is only 
one solution to equation 1151 such that Rt < Rco- A real 
system may deviate slightly from our simple picture (e.g., 
of an unperturbed dipole field), leading to an uncertainty 
in the exact value of Rt. However, due to the steep radial 
dependence of drm relative to dra, the location of Rt should 
not be significantly affected. For the system plotted in Figure 



^ Note that Rt is related to the 'fastness parameter,' uj, in X-ray 
pulsar literature: ui = (Rt/Rco)^^^ ■ 
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log(f) 

Figure 5. Predicted location of Rt/Rt as a function of log(/), 
for the fiducial CTTS system. The dashed lines represent the 
prediction from equation 1151 for /3 = 0.1, 0.5, 1.0, and 2.0, as 
indicated, and all have 7c — » oo so that the magnetic field is 
assumed to remain closed. The solid line is for /3 = 0.1, but with 
7c = 1, so that the field is partially open. The system with 7c = 1 
switches to state 1, for log{/) < — 1.5, and Rt is then predicted 
by eauation ll8l The dotted line shows the location of Rco (eq. |^. 



|1] the solution to eguation llSl is Rt/Rco ~ 0.915, represented 
by the vertical dotted line in the Figure. 

For a system in state 1 (with Rt < Rin), eauation ll5l is 
not valid because the field will open (see t|2.1.31 . A substi- 
tution of Rt < Rin in equation llSl indicates that the system 
will be in state 1 if 

/<(l-/37c)(7c^)-'/', (17) 

and it will be in state 2 for any larger /. Note that condition 
1171 can never be satisfied if Pjc > 1 (since Rin is then unde- 
fined), and so the system would then always exist in state 2. 
For the more likely case that /3yc ^ 1, state 1 is a possible 
configuration of any system. 

To determine -Rt in state 1, instead of using equation 
1121 for the differential magnetic torque, one must consider 
the maximum possible drm in order that the field remains 
closed. This is determined by using 7 = —7c in equation |21 
By setting this equal to equation 1131 one finds 

Rt = {jciPf^'R, = {2-i,f'\GM,y^/\M.Y^'''^i''/\ (18) 

Note that this equation does not depend on the rotation 
rate or radius of the star and has the same dependences 
on other system parameters as in many previous theor eti- 
cal works (e.g., Davidson & Ostrikcr 1973 ; Ghosh & L ambI 
ri979l:IShu et al.lll99l^ . Since Rt < Rin in state 1, the field 
lines will be open inside Rco- Thus, in state 1, the stellar 
field connects only to a small portion of the disc near Rt, 
from where a funnel flow originates, and all exterior field 
lines are open, as shown in the top panel of Figure |3 State 
1 is discussed further in section 

Figure El shows the predicted location of Rt in units of 
Rt for the adopted fiducial parameters, as a function of the 
logarithm of the spin rate /. For reference, the dotted line 



shows the location of Rco (eq. |5J . The dashed lines show Rt 
for systems with 7c cxd and (3 — 0.1, 0.5, 1.0, and 2.0, as 
indicated in the Figure. The solid line shows the /3 = 0.1 
case, but with a more realistic value of 7c = 1. Note that 
Rt is always less than Rco, and that both decrease as / 
increases. For stronger field coupling (smaller /3), the field is 
more strongly twisted, and so Rt is closer to Rco- 

All of the dashed lines in Figure |^ represent systems 
with 7c — > 00, in which the field remains closed for arbitrar- 
ily large magnetic twist. Thus, these systems are always in 
state 2, and the dashed lines are everywhere given by equa- 
tion^] On the other hand, the solid line represents a system 
with 7c = 1, so it is in state 2 only when />0.03 feq. I17II . 
For smaller /, it is in state 1, and Rt is then determined by 
equation 1181 Since equation 1181 is independent of /, state 1 
is represented by the constant value of Rt « 9.17?, in the 
Figure. Real systems will likely have 7c = 1, in which case 
the value of Rt ~ 9.1/?, represents an upper limit for the 
fiducial CTTS system, regardless /3. 

2-3.2 Accretion torque 

We assume that accreted disc material is quickly integrated 
into the structure of the star and the accreted angular mo- 
mentum is redistributed into the stellar rotation profile. So 
there is a torque on the steadily accreting star that is given 
by Ta = Ma.[ldiRt) — L], where ld{Rt) is the specific angular 
momentum of the disc material at Rt and /, is that of of the 
star. Combined with equation and assuming solid body 
rotation of the star, this becomes 

Ta = liUGMtRtf'^ [{Rt/R.f'^ - fc'/] , (19) 

where k is the normalized radius of gyration [k'^ ~ 0.2 for 
a fully convective star; AC96). This formula is valid for a 
system in either state 1 or 2. 

The term in square brackets in equation 1191 is dimen- 
sionless and compares the accreting angular momentum 
(first term) with how much the star already has (second 
term). Note that the first term will always be greater than 
or equal to one, while the second term has a maximum value 
of k^ (when / = 1). Thus, the second term is usually negli- 
gible (particularly when / ^ 1). 

2-3.3 Magnetic torque 

When in state 2 (i.e., when Rt > Rin), the stellar field con- 
nects to a significant portion of the disc, and one can in- 
tegrate equation 1121 over the connected region, from Rt to 
Rout, to obtain the total magnetic torque on the star, 

1 ji^ 

Tm ~ -T-r g [2(i?co/7?out)'^^^ — {Rco / Rout)^ 
op tico 

-2{Rco/Rtf'^ + {Rco/Rtf] - (20) 

This torque is independent of the detailed structure of the 
Keplerian disc. Also, equation 1201 includes the dependence 
of the magnetic torque on the field topolgy via the variable 
Ront- For example, the spin-down torque (found by setting 
Rt ~ Rco) exerted by field lines connected out to -Rout ~ 
2.4-Rco is one half of the spin-down torque for -Rout oo- 
Similarly, for -Rout ~ 7.2-Rco or -Rout ~ 34-Rco, the spin- 
down torque is 90% or 99% (respectively) of the spin-down 
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Figure 6. Same as Figure |1] except tliat tiie equilibrium spin 
period is 3.3 days (iJco ~ 3.7R*), and 7c = 1, so that the magnetic 
field is not connected for r > -Rout ~ 1.6i?co- Also, i?t ~ 0.97i?co- 



Figure 7. Logarithm of the ratio of real to maximal spin-down 
torque ("eg. 1221 . as a function of log{/37c). Dotted lines represent 

X = /37c and x = (/37c)"^- 



torque for -Rout oo. It is evident that, even when -Rout is 
large, most of the spin-down torque comes from field lines 
connected not too far from J?co- This is simply because the 
differential magnetic torque (eq. I12II becomes very weak far 
from the star. Thus the typical assumption of -Rout — > oo 
is not significantly effected by the fact that real discs have 
finite radial extents, so long as they reach to several times 

Above, we have taken -Rout as arbitrary, but our goal in 
this paper is to consider the opening of the field from dif- 
ferential rotation, so -Rout is then given by eauation llOl and 
the preferred formulation of the magnetic torque becomes 



(21) 



"""^ ^ 3^fe [2(1 + /37c)-' ~ (l + /37c)-' 

-2(i?co/i?t)''" + (i?co/i?t)'] . 

This is exactly the solution found by AC96 for the special 
case of /3 = 1 and 7c ^ oo (so that -Rout oo), but our 
formulation includes the effect of field opening via differen- 
tial rotation, in which case, the field topology depends on 
the diffusion parameter /3. The total magnetic torque on the 
star can be either positive or negative, depending on the 
size of the connected region inside 7?co, compared to the 
connected region oustide -Rco- In the next section f2.4l we 
will show that, for reasonable values of /3 and 7c, -Rout is 
very close to 7?co, and the spin-down torque is significantly 
effected. 



2.4 Effect of opened field 

Figurel^shows the differential torques in the fiducial CTTS 
system with /3 = 1, as in Figure |1| However, unlike Figure 
0] here we show the system for 7c = 1, so that the field 
is open beyond -Rout ~ 1.59-Rco- The star now rotates sig- 
nificantly faster, with a period of 3.3 days (/ « 0.14), and 
-Rt « 0.974-Rco- A comparison between Figures 2] and |H] illus- 
trates the effects of varying field topology on the differential 
torques in the star-disc system. 



There are some interesting differences between Figures 
2] and |S| Most notably, the differential magnetic torque in 
Figure|S|abruptly goes to zero at the location of -Rout, due to 
the assumption that there is no torque on the star from the 
disc where the field lines are open. It is evident that a more 
open topology results in a smaller connected region, which 
leads to a net (integrated) spin-down torque that is smaller 
than for the completely closed topology. Thus, a more open 
topology results in a faster equilibrium spin rate, as can be 
seen by comparing the stellar spin rate of 6.0 days for Figure 
|l|with 3.3 days for Figure |S| 

To quantify the effect of a more open topology on the 
magnetic torque and to determine the dependence on /3, 
we first consider only the portion of the magnetic torque 
that acts to spin down the star by setting Kt ~ Rco in 
equation 1211 For normalization, we use the maximum spin- 
down torque, t^I^^ (eq. We define the ratio of the true 
spin-down torque to this maximum torque as x = Tin{Rt ~ 
■Rco)/rt™'L'i, which is given by 



X = (/37c)~' [l + (1 + /37c)-' - 2(1 + /37c)" 



(22) 



and only depends on the parameter /37c. It is immediately 
evident that for /37c = 1, x = 1/4, so the spin-down torque 
is four times less than used by AC96, when one considers a 
more realistic magnetic topology. 

Figure |7| illustrates the dependence of the torque ratio 
X on /37c. It decreases as P'yc'^ for /3 3> 1 and increases 
as /37c for (3 <^ 1 (as revealed by Taylor expansion of eq. 
I22II . The limiting behavior of x is indicated by the two dot- 
ted lines in the Figure. This behavior can be understood as 
a competition between two effects: In the strong magnetic 
coupling limit (/? <C 1), the field topology becomes more 
open for smaller /37c, reducing the spin-down torque. In the 
weak coupling limit (/3 ^ 1), the topology is largely closed, 
but the twisting of the field lines is smaller for larger /3, 
which reduces the differential magnetic torque at all radii 
(drm (X /3~^). These two effects conspire to give a maximal 
value of X for the special case of /37c — 1. For the more 
likely value of /3 = 10~^, the spin-down torque is two or- 



10 Matt & Pudritz 



T — ' — ' — ' — ' — I — ' — ' — ' — ' — r 




log(f) 

Figure 8. Net torque, Ta + Tm (eqs. 1191 and I21i . on the star as 
a function of the logarithm of the fractional spin rate, for the 
fiducial CTTS system. The four solid lines represent (/3,7c) = 
(l,oo), (1, 1), (0.1, 1), and (0.01, 1), as indicated, corresponding 
to different field topologies, ranging from completely closed to 
more open. A system will be in spin equilibrium when the net 
torque equals zero (dotted line). 

ders of magnitude lower than used by AC96. While it may 
at first seem surprising that strong magnetic coupling leads 
to weaker spin-down torques on the star, further reflection 
reveals that this is necessarily true, since stronger coupling 
leads to stronger twisting, which further disconnects the star 
from the disc. 

Finally, equations 1191 and 1211 can be used to calculate 
the net torque on the star. Figure|H|shows this net torque as 
a function of log(/) for the fiducial CTTS parameters and 
for different values of /3 and 7c. For each case, we calcuate 
the net torque as follows: First, for a given value of f3 and 
7c, and for each value of /, we use eauation llTl to determine 
whether the system is in state 1 or 2. We then find the 
location of the truncation radius, using equation 1151 if the 
system is in state 2, or equation 1181 if in state 1. Finally, we 
calculate the integrated torques Ta (eq. I19II and (eq. ITTI 
if in state 2; Tm = 0, if in state 1). As discussed in section 
12. .'-{.II only cases with Pjc < 1 can be in state 1, so only those 
cases show a transition at log / « —1.45 in Figure |5] The 
Figure also shows the effect of the field topology on the net 
torque on the star from the disc. It is evident that, when the 
magnetic field is partially open (7c — 1), the net torque is 
larger than for the case where the field is everywhere closed 
(7c — 00). The spin rate at which the net torque on the 
star is zero indicates the equilibrium spin state, which is the 
topic of the next section. 



3 THE DISC-LOCKED STATE 

The general theory presented in section |5| enables one to 
calculate the net torque (ra + Tm) on the star for any ac- 
creting system with known Af, , i?, , /i. Ma, and Q, (one 
must also adopt a value for 7c and /3). The system is sta- 
ble, in that a positive torque spins the star up, and a faster 



spin reduces the total torque. Conversely, a negative torque 
spins down the star, and the torque increases (becoming less 
negative) for slower spin. Therefore, in a system where the 
other parameters are relatively constant, the spin rate of 
the star naturally adjusts to an equilibrium state in which 
Ta + Tm = 0, kuowu as the 'disc-locked' st ate (e.g., K91; 
ICameron fc Campbelj[l99i IShu et al.lll994l : AC96). Since 
the only torques that spin down the star originate along 
field lines that connect to the disc outside Rco, systems in 
their equilibrium state must be in state 2 {Rt > Rin, see Fig. 
1^1. In this section, we show that both Rt and Q* in the disc- 
locked state are significantly affected by the field topology 
and thus have a strong dependence on the magnetic diffusion 
parameter /3. 

3.1 Truncation radius in the disc-locked state 

The disc-locked state is defined by the condition, Tm ~ — Ta. 
Thus, by combining equations 1191 and 1211 and using equa- 
tions |21 and ^| to eliminate /, this condition can be rear- 
ranged to be 

j^(/37c)-(i^co/i^0cT _^ (23) 
{Rco/Rt)%'[l-iRco/Rt)%'] ' 

where 

;f(/37c) = 2(1 -f /37c)-' - (1 + /37c)"', (24) 

and the subscript 'eq' refers to the value in the disc-locked 
state. In deriving eauation l28l we have ignored the term pro- 
portional to k^f in eauation ll9l as justified in the discussion 
following that equation. The function K{f3'yc) characterizes 
the topology of the field in a sense that, when /37c varies be- 
tween and oo, K varies between 1 (completely open field) 
and (completely closed). Equation 1231 has exactly one so- 
lution such that {Rt/Rco)cq < 1 (which is the only physical 
solution) for any given /37c > 0. 

The location of Rt for accreting sy stems is, in principle , 
an observable parameter. For example. iKenvon et alj lll996l) 
used a magnetic accretion model to predict infrared excesses 
in CTTS's and then to determine the value of Rt/Rco for a 
sample of stars in the Taurus- Auriga molecular cloud. The 
value of {Rt/Rco)cq predicted by eq nation 1231 represents the 
value for a system that is disc locked. 

The solid line in Figure |U] shows {Rt/Rco)cq as a func- 
tion of log(/37c). When a closed field topology is assumed 
(7c — > oo), (/?t/i?co)cq ~ 0.915. However, for the more rea- 
sonable value of 7c — 1, {Rt/Rco)aq increases (approaching 
unity) as the magnetic coupling bec omes st r onger (smaller 
/3) . This confirms the conc lusi on of IWan d il995l and see 
ICameron fc Campbedll993l: lYillTggl 119951) that any disc- 
locked system will have (i?t/i?co)cq ^ 0.9, and we find that 
the effect of a more open field topology is to significantly 
increase this value. 

The Figure also indicates the three possible states of the 
system (discussed further in H4.2II . determined by the loca- 
tion of Rt relative to Rm and Rco (dashed lines). A system 
that is disc-locked is always in state 2. If a system is observed 
with Rt /Rco larger than the solid line in the Figure, the star 
should be spinning down. Conversely, if Rt/Rco is smaller 
than the solid line in the Figure, the net torque from accre- 
tion and from field lines connecting the star and disc will 
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Figure 9. Location of Rt/Rco in the disc-locked state (solid 
line; eg. I2ijl . as function of log(/37c). The dashed lines represent 
Rin/Rco and the location where Rt/Rco = 1- Any system with 
Rt > ^?co is in state 3 (dark grey shaded region) , while any system 
with _Rt < Rin 

is in state 1 (light grey). Purthermore, if Ri/Rcq 
is anywhere below (above) the solid line, the net torque from the 
disc will spin the star up (down). 



Figure 10. Spin rate factor C(/3,7c) (eg. 1261 as a function of 
log(/3). The dotted line is C(/3) for 7c — » oo and represents the 
assumption of a completely closed magnetic topology. The solid 
line shows the more realistic case with a partially open topology 
{■jc = !)■ The dashed lines represent s C ~ 1.15 from K91. The 
Figure is from lMatt fc Pudritd i2004) . 



act t o spin the star up. It is interesting that iKenvon et alJ 
(|l996) found typical values of Rt/Rco in the range of 0.6 to 
0.8 for the stars in their sample. If true, these stars cannot 
be in spin equilibrium, unless they feel significant spin-down 
torques other than those from field lines connecting them to 
their discs. Furthermore, if /3yc < 0.3 is appropriate, the 
stars in their sample should exist in state 1. 

3.2 Stellar spin rate in equilibrium 

Now that we can calculate {Rt/Rco)oq via equation 1231 we 
rewrite equation 1151 to find the fractional spin rate of the 
star in equilibrium 

/cq = C(/3,7c)(2/V')'/', (25) 
where 




Since {Rco/Rt)cq depends only on f3^c (via egs. 1231 and I24II . 
the dimensionless function C(/3, 7c) depends only on /3 and 
7c. We can also combine equations and l25l to find the 
equilibrium angular spin rate of the star 

QT = C{f3, 7c)Ma''^ (GM,)^/^ /i-^/^ (27) 

This equation has the same dependence of ^1%'^ o n Ma, M,, 
and u as e quation 3 of K91, and a s in the theory of lShti et alJ 
Jl994l and lOstriker fc S hu 1995). The only difference is the 
value of the f actor C used in the var ious theories. K91 used 
C ^ 1.15, and lOstriker fc Shd(|l995^ '' found C « 1.13. How- 
ever, our formulation of the problem allows us to determine 



the effect of the field topology on the equilibrium spin rate, 
via the function C(/3, 7c), for arbitrary values of the diffusion 
parameter /3. 

Figure \Wl reveals the dependence of C(/3, 7c) on /3 for 
the two values of 7c we have considered throughout. For 
7c ^ cxD, C(/3) ^ 1.59/3^^^, which is represented by the dot- 
ted line in the Figure. The solid line shows C(/3, 7c) for the 
more realistic value of 7c = 1 and illustrates the effect of 
a reduced magnetic connection to the disc. For comparison, 
the dashed line shows the spin rate factor used by K91, 
which also roughly represents the typical factors of order 
unity used in most star-disc interaction models. 

The dotted line in Figure 1101 represents the assump- 
tion that the magnetic field everywhere connects to the 
disc, regardless of the field twist. In that case, the mag- 
netic torque increases with with decreasing /3, since the field 
then becomes highly twisted, and so the prediction is that 
SI*'' oc Z?"'^^. However, when one considers that dipole field 
lines will become open when largely twisted, the torque has 
a maximum value for /3 = 1 and decreases for any other (3 
(see il2.41 . Correspondingly, the solid line in Figure 1101 has 
a minimum value for /3 = 1. This minimum value represents 
the 'best case' for disc locking, and even at this location, C 
is a factor of 1.8 larger than the value for 7c ^ oo and 2.5 
times larger than used by K91. For the more likely value of 
/3 = lO"'^, the predicted equilibrium spin rate of the star 
is more than an order of magnitude faster than prediced 
by any other model. Note that the spin rates of the system 
plotted in Figures |1] and |S| were chosen to be in equilibrium, 
and a comparison between the two Figures shows the effect 
of fi eld topology for the = 1 'best case.' 

iMatt fc Pudritd ll2004l) applied the analysis presented 



^ While it is interesting that our eguation 1271 resembles the for- 
mulation of lOstriker &: Stnl ^^9^, their assumed magnetic field 



geometry is different than ours, so the comparison of C values 
should not be taken too seriously. 
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Figure 11. Equilibrium spin period as a function of observable parameters Ma (upper left), (upper right), M, (lower left), and i?* 
(lower right). In each panel, all parameters are held fixed at the fiducial value except for that which is plotted along the abscissa. The 
vertical dotted line in each panel shows the fiducial value that is held fixed in the other panels. The solid line in each panel represents 
the assumption that the magnetic field topology is completely closed (7c 00) and /3 = 1. All of the broken lines represent a partially 
open topology (7c = 1) for /3 = 1 (dash-dot lines), 0.1 (dashed lines), and the more likely value of 0.01 (dotted lines). 



in this section to the CTTS BP Tau, which is one of the 
few stars for which all of the relevant system parameters are 
known (or well-constrained) . They argued that the existence 
of slowly rotating, accreting stars, such as BP Tau, cannot be 
explained by a disc-locking scenario. To further illustrate the 
effect of field topology on the predictions of disc-locking, and 
to apply the analysis to all CTTS's, we have plotted Figure 
111! The Figure shows the predicted spin period for a wide 
range of observable parameters. The spin period is given by 
27r/r2°'', and we have used the relationship /i — BtR^. 

The solid lines are for 7c — > 00 and j3 — 1, and so they 
represent the 'standard prediction' by previous models. The 
broken lines take into account that some of the field should 
be open (7c = 1) for the three different values of /3 = 1, 
0.1, and 0.01. Note that [3 = 0.1 predicts the same period as 
for /3 — 10 (due to the approximate symmetry of the solid 
line in Fig. llOH . and /3 — 0.01 corresponds with /3 = 100. 
It is evident that, even in the 'best case' (/3 = 1) for the 
disc-locking scenario, the effect of a more open topology is 
to reduce the equilibrium spin period by a factor of two, 
compared to the closed field assumption. Given the uncer- 



tainties in some of the observed parameters, it may not yet 
be possible to constrain the predicted period to within a fac- 
tor of two. In particular, a difference of a factor of two in the 
predicted spin period could result from an error of a factor 
of 5.0, 2.2, 2.6, or 1.3 in the observed parameters Afa, B,, 
A/*, or R,, respectively. However, for the more likely value 
of /3 = 10~^ (dotted lines in Fig. IllH . the predicted spin 
period is an order of magnitude lower than the 'standard 
prediction,' which cannot be reconciled by observational er- 
rors. 

3.3 Time to reach equilibrium 

It is important to determine how quickly the star will spin 
up or down to reach the equilibrium state. Rather than fully 
solving the time-dependent problem, which should also in- 
clude thesmn-up due to the contraction of the protostar 
(e.g.. lYll994l) . one typically estimates a characteristic spin- 
down time using the angular momentum of the star, L,, 
divided by the net torque on the star. To be more precise, 
and for arbitrary spin rates, one should replace L* with the 
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Figure 12. Time to reach equilibrium, ispin feQ. I28i . as a function 
of the spin rate fraction /, for the fiducial CTTS system. The four 
solid lines represent (/3,7c) = (1, oo), (1, 1), (0.1, 1), and (0.01, 1), 
as indicated. The vertical dotted lines indicate the equilibrium 
spin rate, for each case. 

difference between the current L* and the value of L* for 
the equilibrium spin rate. Assuming solid body rotation of 
the star, the characteristic time to reach spin equilibrium is 
then 

\ Ta + Tni J 

which corresponds to a spin-up (down) time for a star cur- 
rently spinning slower (faster) than Q!^^. If ts-pi-a is long 
compared to the time-scale for other system parameters to 
change (e.g., compared to the lifetime of the disc), the star 
is unlikely to ever be in a spin equilibrium state. 

Figure WI\ shows ts-pi-n for the adopted fiducial parame- 
ters, as a function of the spin rate fraction / and for differ- 
ent values of /3 and 7c. It is evident from the Figure that 
tapin generally decreases for increasing /, since faster spin 
means i?co is closer to the star so the magnetic torque will 
be stronger (whether it spins the star up or down) . There is 
an exception to this when the star spins much slower than 
/cq. For example, for the case with (/3,7c) = (0.01, 1), when 
/ < 0.2, tapin decreases with decreasing /. This can be under- 
stood, since then Rt/Rco decreases rapidly with decreasing 
/ (see Fig. ISJ, leading to much stronger magnetic spin- up 
torques. Note also that the cases with (/3, 7c) = (0.01, 1) and 
(0.1, 1) are in state 1 for /< 0.036 and 0.033, respectively, 
in which Tm ~ 0- 

For the 'standard' prediction with (/3, 7c) — (1, <x), 4pin 
is always less than 5 x 10* years, which is much shorter 
than the expected di sc lifetime of more than 10® years 
jMuzerolle et al.ll2000l) . However, when one considers the 
effect of a more open field topology (7c ~ 1), the mag- 
netic torques is reduced, and so tapin is longer and increases 
when /3 decreases. For the likely case with (/3, 7c) = (0.01, 1), 
tspin ~ 4 X 10^ years. This is still relatively short compared 
to the expected disc lifetime. Therefore, we agree with pre- 
vious authors (e.g., K91; AC96) that systems such as those 
considered above should exist near their equilibrium spin 



states throughout most of their accretion lifetimes, but only 
if /3>0.01. However, the effect of a more open field topol- 
ogy is that the equilibrium spin rate is much faster than 
previously predicted. 

We note in passing that the characteristic time fspin we 
have calculate d here is s i gnifica ntly shorter than recently 
calculated by iHartmaiml i2002h . Hartmann's estimate as- 
sumes that the upper limit to the spin-down torque on the 
star is equal to A/a(GM,iiout)^'^^. However, as discussed in 
section 031 this is the torque necessary to provide a steady 
accretion rate of Ma through the radius Rout- In order for 
the disc to exert a spin-down torque on the star (via the 
magnetic connection), it must provide a torque in addi- 
tion to Ma,{GMtRo\itY^^ , requiring the disc to have a dif- 
ferent structure than in the absence of a stellar field (see 
iRappaport et af]|2004^ . It is not yet clear to what extent 
the disc can be restructured (before accretion will cease), 
and so Hartmann's estimate does not represent a true limit. 



4 DISCUSSION AND CONCLUSIONS 

We extended the standard picture of the interaction of mag- 
netized star with a steady-state accretion disc. Our more 
comprehensive formulation of this problem allows us to de- 
termine the location of the disc truncation radius Rt and 
calculate the torque on the star for a system with arbitrary 
values of Afa, M*, i?*, B, , and /3, which parametrizes 
the coupling of the magnetic field to the disc. We consider 
only two sources for the torques: (a) torque from the angular 
momentum deposited by accretion of disc material from Rt , 
and (b) torques exerted by field lines connecting the star to 
the disc over the region from _Rt to -Rout. 

Our model resembles several previous studies (e.g., 
AC96), except that we have now determined the dependence 
of the torques on the magnetic coupling to the disc. Specif- 
ically, the difi'erential rotation between the star and disc re- 
sults in a largely open topology (e.g., UKL), so the size of the 
region of the disc that is magnetically connected to the star 
is smaller (i.e., -Rout is smaller). Thus, when one considers 
this effect, the magnetic spin-down torque on the star is less 
than if one assumes the field remains everywhere closed. The 
strongest spin-down torque occurs for intermediate magnetic 
field coupling {(3 = 1), in which case the spin-down torque 
is a factor of four less than for the closed field assumption. 
For strong magnetic coupling (/3 ^ 1), as expected near the 
inner edge of an accretion disc, -Rout is very close to -Rco, 
and the spin-down torque then is proportional to jS. For the 
likely value oi 13 = 0.01, the spin-down torque is 100 times 
less than for the closed field assumption! Furthermore, the 
possibility that field lines may open inside -Rco characterizes 
a new mode (state 1) in the system, in which the star feels 
no spin-down torques from field lines connected to the disc. 
Three possible system states are summarized in section 

We also considered the disc-locked state of the system, 
in which the net torque on the star is zero. A more open 
field topology leads to an equilibrium state that has a higher 
stellar spin rate. The time for a given system to reach spin 
equilibrium is also longer when the field is more open. These 
results apply to any system in which accretion occurs onto a 
magnetized central object. In the general case, not all of the 
system parameters are observationally known. Thus, one of- 
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ten assumes that a particular system is disc-locked, and then 
'tunes' the unknown system parameter(s) to satisfy equa- 
tion |53 We found that, equation 1271 contains the function 
C{P, 7c), which is plotted in Figure lTUl (solid hue), as a func- 
tion of l3. It is evident, that when the coupling of the field 
to the disc is strong (/3 ^ 1) or weak (/3 1), the 'tunable' 
system parameters will be significantly different than for the 
usual assumption that C(/3, 7c) is a constant near unity. In 
particular, systems that are thought to be disc- locked will 
require a larger ^, larger Ql'^, smaller M*, or smaller Ma 
than calculated using previous formulations of equation 1271 
In the next section, we discuss the implications of these re- 
sults, and recent results from the literature, for disc-locking 
in CTTS's. 



4.1 Problems with disc locking for CTTS's 



still controversial ( 


in particualar, see IStassun et al. 


1 119991: 


iHerhstet al .11200(1: 


Istassun et alJl200ll: iHerbst et al.l 


l2QQd), 



so we have taken another look at the problem from a the- 
oretical standpoint. We find that spin-down torques on a 
CTTS are significantly reduced for strong coupling of the 
stellar field to the disc (small /3), resulting in equilibrium 
spin periods as low as a few days or less, for a wide range of 
system parameters (see Fig. llljl . Small values of /3 are likely 
for CTTS's (see >I2.1I| . although /3 is a highly uncertain pa- 
rameter to calculate from first principles and may even vary 
from system to system. Also, as discussed in section [2.1.21 
there may exist special circumstances that allow the field to 
remain connected, but whether these circumstances can ex- 
ist in CTTS systems remains to be shown. Therefore, while 
our analysis of the problem (in does not completely rule 
out the possibility for disc locking in all systems (particu- 
larly for fast rotators) , it significantly reduces the likelyhood 
that disc-locking can expla in the existence of a ccreting stars 
spinning at ~ 10% (e.g., iBouvier et alJll993h of break-up 
speed. There are several other issues in recent literature that, 
when combined with our analysis, cast additional doubt on 
the applicability of disc-locking to the slow rotators ^. 

The most notable observational challenge to the disc- 
locking model is the appare nt lack of str ong dipole fields on 
CTTS's (first suggested bv lSafieil[l993) . It is generally ac- 
cepted that CTTS's have field strengths of a few kiloguass. 
This was predicted by K91, and subsequent observat ions 
iBasri et al.. 1992: Guenthcr et al 1999: Johns-Krull efai] 
I1999I: I.Tohns-Krull fc Valentill200nl: rTohns-Krull et al.ll200lll 
have indeed found a mean field on the surface of the cen- 
tral stars of typically 2 kG. Thus far in our analysis, we 
have adopted a field strength of 2 kC to represent a typi- 
cal CTTS system and guide our discussion. However, strin- 
gent measurements of the mean line of si ght field have been 
carrie d out for t hree CTTS's, BP Tau Oohns-KruU et alJ 
[l99i), TW Hya dJohns-KruU fc Valentill200 J) , and T Tau 



^ Due to the the unique field geometry of the X-wind, in which all 
the stellar field lines are squeezed into the inner edge of the disc 
iShu et al.ll994tl , that model avoids the problem of field opening 
due to differential twisting, as considered in this paper. However, 
the issues discussed after the first paragraph in this section apply 
to all disc-locking models, including the X-wind. 



iSmirnov et al.ir2003l |20o3), and all measurements give an 
upper limit of roughly 200 G for the strength of the dipole 
component of the stellar magnetic field. The measured mean 
fields of 2 kC thus represent a field that is disordered or 
characterized by multip oles of higher order than a dipole 
jjohns-Krull et alll999^ . Such high order fields, even if very 
strong on the stellar surface, decrease in strength too quickly 
with increasing radius to exert significant spin-down torques, 
especially for slow rotators in which Rco is at several stel- 
lar radii. Furthermore, a 200 G dipole field cannot exert 
a significant spin-down torque on a CTTS, even if the field 
connects to the disc everywhere outside J?co. This is evident, 
for example, in the upper right panel of Figure ITTI which in- 
dicates that the equilibrium spin period for a star with such 
a field is less than one day. For the cases with 7c = 1 and 
/3<0.1, there is no equilibrium possible, since the magnetic 
spin-down torques are not strong enough to counteract the 
angular momentum added by accretion, even for maximal 
stellar spin (/ = 1). Also, the time to reach equilibrium 
(for cases in which equilibrium is possible; as discussed in 
increases by an order of magnitude when B, — 200 G, 
compared to 2 kC. 

Second is the issue pointed out by IWand il995ft and 
discussed in section HH. II that stars in their disc- locked state 
must have Rt/Rco^O.9, for a wide range of possible as- 
sumptions in the model. Interestingly, iKenvon et al.l dlQOdl 
concluded that, for their CTTS sample, the typical value 
of Rt/Rco was in the range of 0.6 to 0.8, well below the 
disc-locked value. If true, these stars cannot be disc-locked, 
since the sum of the accretion torque and the torque carried 
by field lines connected to the disc will be positive (spin- 
ning the stars up; see il2.3.3l and 13.111 . Thus, these stars can 
only be in spin equilibrium if they feel significant spin-down 
torques other than those from field lines connecting them to 
their discs. This conclusion does not depend on whether or 
not the stellar field can open, since the calculation of Rt in 
equilibrium also does not. 

F inally, CTTS's may drive stellar winds iSafieij 
'199^, and outflows from the disc are known to ex- 
plain several aspects of observed protostellar outflows (e.g. 
iKonigl fc Pudritdl2000h . Ionized winds escape from regions 
with open magnetic field lines, or they can the mselves open 
the field (e.g., as in th e solar wind: [ Par kerll9'53^ . disconnect- 
ing the star and disc. ISafieil figOS) concluded that CTTS's 
winds should open all stellar field lines beyond roughly 3 
R,. Furthermore, a recent m easurement of rotatio n in the 
jet from the CT TS DG Tau (iBacciotti et al.ll200l and see 
iTesti et al.ll200^ suggests that the low velocity component 
(~ 70 km s~^) originates in the disc from as close as 0.3 AU 
from the star jAnder son et al.ll2003li . This is likely an upper 
limit (|Pesenti et alJl2004^ . and the more ti ghtly collimated , 
high velocity component 220 km s~^: IPvo et al.ll2'003h 
must then originate fr om well within this radius in the disc 
([^i^grggi^^^Jj222^. Theoretical disc- wind models (e.g., 
lK5ni^^^Pu^^3200(l) predict jet speeds of the order of the 
Keplerian velocity from where the wind is launched, and ob- 
served protostellar jets typically travel with speeds of a few 
hundred km s~^ (e.g.. iR.eiourth fc B allv 2 001iV Considering 
a star with M« — 0.5Mq and i?« = 2Rq, and assuming 
Wk = 100 km s~^ at the launch point, this requires the disc 
winds to originate from near GM,/vl ^ 4.8i?.. These ion- 
ized disc winds prevent the star from being magnetically 
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connected to the disc beyond the innermost location of the 
wind launching point (effectively giving an upper limit to 
-Rout), and these winds carry angular momentum from the 
disc, not from the star. As calculated in this paper, a re- 
duced size of the connected region leads to a reduced spin- 
down torque on the star, regardless of the cause of the field 
opening. Protostellar outflows thus provide an independent 
and stringent constraint on disc-locking models. 

We conclude that spin-down torques exerted by field 
lines connecting the star to the disc outside -Rco are likely to 
be much weaker than usually assumed. Therefore, the exis- 
tence of slowly rotating (/ < 0.1, or perhaps higher) CTTS's 
probably cannot be explained by a disc-locking scenario. Ei- 
ther these stars are all in the process of spinning up, or the 
stars feel torques other than those related to a magnetic 
connnection to the accretion disc. Given that the typical 
spin- up times for these systems are short (see tl3.31 , the lat- 
ter possibility appears the most likely. 

4.2 Three states of the system 

We identified three possible configurations of the system, 
determined by the location of Ri, relative to the two key 
radii -Rin fea. llH and Rco- There are two accreting configu- 
rations, which we call states 1 and 2, and one non-accreting 
configuration, state 3. Here, we summarize the conditions 
that determine and characterize each state. 

Figure illustrates the basic magnetic configuration of 
the three states, which may even represent an evolutionary 
sequence for a system with (e.g.) an evolving Ma (see 
Figure |n| shows the location of Rin/Rco (curved, dashed line) 
for various values of /37c, and the horizontal dashed line 
represents the location of Rco- For a system with a given 
value of /37c and for which Rt is determined, the Figure 
indicates which state the system will be in and whether the 
star should be spinning up or down. 

4-2-1 State 1: Rt< Rin (or R'out < Rco) 

This state is a direct consequence of the opening of field 
lines via the differential rotation between the star and disc. 
When Rt is sufficiently less than Rco (i.e., when Rt < Rin), 
the magnetic field becomes highly twisted there, opening the 
field inside Rco and resulting in a highly open field topology. 
Note that state 1, in this context, is only possible for /37c < 
1, since otherwise _Rin is undefined (see Fig.|3]and discussion 
in ilTOl . 

As illustrated in the top panel of Figure|21and discussed 
in section r2. 3. II the stellar field in state 1 connects only to 
a very small, innermost region of the disc near Rt, and all 
exterior field lines are open. The size of the small connected 
region is likely to be determined by dissipative processes 
within the disc, but we do not attempt to calculate this 
here. The location of -Rt is determined by eauation llSI (and 
see Fig. I^J, and accretion onto the star occurs from there. 

A star in this state will always feel a net positive torque 
from the disc, since no field connects outside Rco- Specifi- 
cally, the star is spun up by the accretion torque (ra; eq. 
I19II . and equation 1211 for the magnetic torque is not appli- 
cable. In fact, since Rt increases with stellar field strength, 
and since Ta oc -R^^, the presence of a stellar field actually 



increases the spin up torque on the star, relative to non- 
magnetic accretion. In any case, a system in state 1 cannot 
be in spin equilibrium, unless it receives torques other than 
those considered in this paper. 

Is state 1 a likely, or even common, configuration for 
accreting systems? Figure |S| indicates that, when the mag- 
netic coupling to the disc is strong, the range of possible 
values of Rt for which a system can be in state 2 is signifi- 
cantly reduced. For example, if /37c = 0.01, any system with 
-Rt < 0.993-Rco should be in state 1. The specific conditions 
under which any given system will be in state 1 is given 
by eauation ll7l For illustrative purposes, we can solve this 
equation (and using eg. 1161 for the mass accretion rate. As- 
suming 7c = 1, /3 — 0.01, and / — 0.1, and using the fiducial 
CTTS values (see discussion above >I2.H . one finds that a 
system should be in state 1 if 




This threshold value of is an order of magnitude larger 
than the fiducual value, suggesting that CTTS slow rotators 
will most commonly exist in state 2. 

However, eauation l29l assximes a magnetic field strength 
of 2 kG, and as discussed in section ITTI the stars likely have 
surface dipole field strengths of less than 200 G. This con- 
sideration decreases the threshold value of -Ma by at least 
two orders of magnitude, suggesting that typical CTTS sys- 
tems may exist in state 1. We can also look at this from 
the standpoint of stellar spin, using equation 1171 which in- 
dicates that a system with the adopted fiducial parameters 
(but with B, = 200 G) will be in state 1 if it spins more 
slowly than 26% of break-up speed. (As discussed in section 
12.3.11 this also corresponds an upper limit of -Rt <2.4-R».) 
Thus, if the dipole fields are indeed weak, it is more likely 
that slow rotators, and even some fast rotators, will be in 
state 1. Furthermore, the conclusion ofJCcnvon ct al. ( 199^), 
that Rt/Rco typically ranges from 0.6 to 0.8, suggests that 
CTTS's will be in state 1, as long as /37c < 0.3 (see ^TTl . 

We have thus far considered the opening of field lines 
via the differential rotation, so the existence of state 1 re- 
quires that /37c is significantly less than unity. Given the 
large uncertainty in the value of /3 in real systems, it is still 
not clear whether or not state 1 should be common. From 
the standpoint of torques on the star, the most important 
feature of state 1 is that the stellar field never connects to the 
disc outside Rco- Thus, for the following discussion, we will 
generalize the definition of state 1 to include any magnetic 
configuration in which the star does not connect outside -Rco- 

State 1 is characterized by a large amount of open stel- 
lar field, so it is natural to consider the effects of a stellar 
wind in the magnetically open region. As discussed in sec- 
tion l4.1l a wind can even be responsible for opening the field 
(which does not depend on our parameters /3 and 7c). Thus, 
if a wind (or any other process) keeps the stellar field open 
beyond some radius R'out, and if -Rnnt < Rco, the system will 
be in state 1. For example. ISafieJ (Il998l) concluded that stel- 
lar winds from CTTS's could result in .Rout ^3^*- K true, 
this means that any system rotating more slowly than 19% 
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of break-up speed will have -Rout < Rco (eq. and be in 
state 1. 

There is empirical evidence for systems in state 1 
from some numerical simulations of the star-disc inter- 
action, which usually r epresent systems wi t h <^ 1. 
In the sim ulations of iGoodso n fc Wingled (Il999l) and 
Ivon Rckowsk i fc Brandenbur^ ~ il20o3l . as an example, af- 
ter the initial state, the stellar field never connects to the 
disc outside Rco, even immediately following reconnection 
events. These authors report that the only significant spin- 
down torques on the star come from the open field regions 
(though a stellar wind was not properly included), rather 
than al ong field lines connectin g the stars to their discs (but 
also see iRomanova et "^11200211 . 

It seems that state 1 is a likely configuration for accret- 
ing stars, particularly among slow rotators. Since, in this 
state, the net torque from the interaction with the accre- 
tion disc only acts to spin up the star, stars with long-lived 
accretion phases must somehow rid themselves of this ex- 
cess angular momentum. Stellar winds can exert spin-down 
torques on the star,and if these torques are significant (e.g.. 
Fron t & Pringlc 199^ , the equilibrium spin rate may be sim- 
ply determined by a balance between this torque and Ta. In 
this situation, state 1 could actually represent the expected 
configuration for accreting systems in spin equilibrium. 

4.2.2 State 2: Rin < Rt < Rco 

In this state, the stellar field connects to a finite region of 
the disc between Rt and -Rout, as illustrated in the middle 
panel of Figure |21 This represents the typical configuration 
in many star-disc interaction models, except that the deter- 
mination of -Rout varies between models. The location of -Rt 
is determined by equation^] (and see Fig.|5Jl, and accretion 
onto the star occurs from there. 

The star is spun up by the accretion torque (eq. I19II 
and magnetic torques (eq. VZll from field lines connected to 
the region of the disc between -Rt and -Rco and spun down 
by magnetic torques from field lines connected between -Rco 
and -Rout. Therefore, a system can exist in an equilibrium, 
disc-locked state, in state 2, in which the net torque on the 
star is zero, and the spin rate of the star then correlates with 
accretion parameters (see 

When one considers that the differential rotation de- 
termines -Rout (via eq. IIUL both fll'^ and {Rt/Rco)cq are 
larger for smaller /3. Also, as shown in Figure |5] the range 
of (non-equilibrium) values of -Rt that exist in state 2 be- 
comes narrower as f3 decreases. For the strong coupling 
case of /37c — 0.01, a system can only be in state 2 if 
0.993 < Rt/Rco < 1.0. So for strong coupling, it is unlikely 
that a given system will exist in state 2, unless it is very 
near its disc-locked state, which then requires a fast stellar 
spin. 

Another intriguing effect of a largely open field topology 
is that, in order for a system to be disc-locked, the differen- 
tial magnetic torque in the disc drm must be stronger (com- 
pared to the completely closed assumption) in order to make 
up for the decreased size of the connected region (compare 
Figs. 2] and |SJ . When the connected region is very small 
(i.e., for small /3), drm is very large, which ought to have 
a significant effect on the disc structure there. There may 
be a physical limit, beyond which the disc cannot respond. 



and accretion will cease (see discussion of state 3, below). 
This could possibly l ead to a time-dependent process (e.g., 
ISprui t fc Taam' 1993), perhaps analogous to the simulations 
of (e.g.) Goodso n ct al.. (.1999. though their stellar field does 
not connect outside -Rco), and in which there may still be a 
time-averaged net torque of zero. 

4.2.3 State 3: Rt > Rco 

Since no accretion onto the star occurs, state 3 may char- 
acterize the non-accreting, weak-line T Tauri phase of pre- 
main sequence stellar evolution. Also, since the disc does 
not extend inside -Rco, the star feels no spin-up torques, only 
spin-down torques, and so it cannot exist in spin equilibrium 
(Sunvacv fc Shakura 1977b). A possible magnetic field con- 
figuration of this state is illustrated in the bottom panel of 
Figure 

It is not yet clear under which conditions a system 
will be in this state, though the relative values of differ- 
ential torques (or stresses) in the disc are certainly impor - 
tant. Some authors fe.g.. IWanj|l995l: IClarke et alJll995^ 
have speculated th at state 3 occurs when I drm/d Tal becomes 
greater than one iCameron fc CampbelJ [l993 suggested a 
value of 2) anywhere outside -Rco. In this work, we have as- 
sumed that the disc will be structured such that e quation 
ITU is satisfied (see i|2.3l and iRappaport et alJliool) . How- 
ever, this assumption must eventually break down for large 
enough /, large enough fj,, or small enough -Ma. 

The general question of what conditions govern a sys- 
tem in state 3, to our knowledge, remains an unanswered 
astrophysical problem. It is not clear what will determine 
the location of -Rt in state 3 (since neither of equations 
1151 and 1181 is then valid) . In addition to magnetic torques, 
outflows and /or radiation from the star may be important 
Jjohnstonell995l) . Understanding this state is probably rele- 
vant to understanding the transition from classical to weak- 
line T Tauri phases, and it may even have furt her impli- 
cations for gas giant plan et formation/migration iLin et alJ 
liggdiTrilling et all2002ll and for the ultimate dissipation of 
the gas disc. 

4.3 Conclusions 

We have considered that the opening of magnetic field lines 
expected from differential rotation in the star-disc interac- 
tion results in a largely open field topology. This significantly 
alters the torque that a star receives from its accretion disc, 
compared to previous models that assume a closed field. Our 
main conclusions from this work are the following: 

(i) This more open field topology resuts in a weaker spin- 
down torque felt by the star from the disc (' i|2.4ll . The 
strongest possible torque occurs for intermediate magnetic 
coupling to the disc. Stronger coupling, as expected near the 
inner edge of the disc, results in a spin-down torque that is 
more than an order of magnitude below the torque found 
for the closed field assumption. 

(ii) In the disc-locked, spin equilibrium state, this results 
in a stellar spin rate that is much faster than predicted by 
previous models f il3.21 . 

(iii) We have identified and discussed three possible mag- 
netic field configurations in magnetic star-disc systems 
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( i|4.2t . The three configurations could represent, for exam- 
ple, an evolutionary sequence for a system with a gradually 
decreasing mass accretion rate (or, e.g., a gradually increas- 
ing stellar spin rate). Our conclusions about each state, in 
the context of T Tauri stars, are the following: 

(a) Assuming strong magnetic coupling to the disc, 
slowly rotating CTTS's should be in state 1 if Ma ^ 5 x 
1O"^M0 yr-i (eq.HHfor = 200 G). Since typical ac- 
cretion rates are higher than this, state 1 may represent a 
common configuration in these systems. In this state, the 
star feels no spin-down torques from the disc. However, 
if spin-down torques from (e.g.) a stellar wind are signifi- 
cant, stars may be in spin equilibrium in state 1, though 
they should not then be considered 'disc locked.' 

(b) State 2 is the typical configuration assumed in star- 
disc interaction models. For strong magnetic coupling in 
the disc, or if stellar or disc winds are important, we find 
that a system can only be in state 2 under special circum- 
stances. In particular, the accretion rate must be lower 
than for state 1. 

(c) State 3 likely represents the non-accreting, weak- 
line T Tauri phase. Given that the accretion disc can 
restructure itself in response to (e.g.) external magnetic 
torques, it it not yet clear when a system will transition 
into this state. This important evolutionary phase requires 
more theoretical study. 

(iv) These considerations, and additional issues from the 
literature, suggest that slowly rotating CTTS's probably 
cannot be explained by a disc-locking scenario ( t|4.H . 

(v) If slowly rotating CTTS's are in spin equilibrium, 
then another spin down torque must be active in the system. 
We suggest that this might arise from magnetized stellar 
winds. 
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